String Breaking in Non-Abelian Gauge Theories with Fundamental Matter Fields 



Owe Philipsen 1 and Hartmut Wittig 2 * 
1 Institut fur Theoretische Physik, Philosophenweg 16, D-69120 Heidelberg, Germany, 
email: o.philipsen@thphys.uni-heidelberg.de 
2 Theoretical Physics, 1 Keble Road, Oxford 0X1 3NP, UK, 
email: h. wittig 1 ©physics, oxford, ac. uk 



o\ ■ 
o\ : 

On . 

3 ■ 

(N ' 
> 

O . 
(N ■ 

O ■ 

: 

o . 
oo ■ 

5 : 
> . 

I 



X 



We present clear numerical evidence for string breaking 
in three-dimensional SU(2) gauge theory with fundamental 
bosonic matter through a mixing analysis between Wilson 
loops and meson operators representing bound states of a 
static source and a dynamical scalar. The breaking scale is 
calculated in the continuum limit. In units of the lightest 
glueball we find r\, mc ~ 13.6. The implications of our re- 
sults for QCD are discussed. 

PACS: ll.15.Ha, 12.38.Aw, ll.10.Kk 



It has been a long-standing problem in QCD to detect 
"string-breaking" , i.e. the breakdown of linear confine- 
ment, by means of computer simulations. String breaking 
is expected to occur when the energy of a gauge string 
between static colour sources separated by a distance r 
becomes as large as the energy required to produce a 
light quark-antiquark pair, each of which is subsequently 
bound to the static colour sources. Thus, for distances r 
below a typical scale n, the static quarks support a gauge 
string, resulting in a linearly rising potential, whereas for 
r > rb the ground state potential is bounded by the en- 
ergy of two static-light "mesons", V(r) ~ Emm- The 
usual computational strategy is to extract V(r) from Wil- 
son loops calculated in Monte Carlo simulations and to 
look for signs that the linear rise in V(r) turns into a con- 
stant behaviour as the separation r is increased beyond 
the breaking scale. Despite recent efforts hard ev- 

idence for string breaking has not been presented so far. 

It has been suggested that the apparent failure to de- 
tect string breaking is due to insufficient overlap of Wil- 
son loops with the two-meson state 0. Indeed, a recent 
analysis of the problem using strong coupling ideas ||] 
suggests that string breaking is a mixing phenomenon, 
involving both the string and the two-meson state. Thus, 
in order to confirm the mixing picture, the conventionally 
used Wilson loops have to be supplemented by explicit 
two-meson operators. However, given the complexity and 
computational cost of full QCD simulations, it will take 
a big effort to obtain an unambiguous signal for string 
breaking, even if mixing is taken into account. 

In this letter we study the mixing between the string 
and a two- "meson" state in a simplified theory, the three- 



dimensional SU(2) Higgs model. We present clear evi- 
dence for string breaking from Monte Carlo simulations 
and essentially confirm the mixing scenario ofref. @. Al- 
though the three-dimensional SU(2) Higgs model is usu- 
ally applied to physical contexts other than the strong 
interactions, its phase diagram has a confinement region 
that shares a number of features with QCD, such as a 
linearly rising potential of static colour charges and its 
eventual screening through matter pair creation. As in 
QCD, the screening phenomenon has so far not been ob- 
served in simulations employing Wilson loops only |J . In 
our previous studies of the model ||-[l0| we have investi- 
gated a possible mechanism how the confining properties 
of the model are lost in the Higgs region through flux 
tube decay Furthermore, an attempt to estimate the 
screening length by using non-local gauge-invariant op- 
erators has been described in p3| . 

Bosonic matter fields imply an enormous simplifica- 
tion of the computational effort, whilst preserving the 
underlying mechanism for string breaking to occur. The 
presence of a weak or even moderately strong scalar self- 
coupling, which has no analogue in QCD, is not expected 
to change the dynamics of gauge and scalar fields so dra- 
matically that the interpretation of results in a QCD con- 
text is no longer possible . 

The study of the model in three dimensions also offers 
a number of advantages: first one can simulate relatively 
large volumes, so that large separations of colour charges 
can be probed efficiently without suffering from signifi- 
cant finite- volume effects. Second, due to superrenormal- 
isability, the curves of constant physics can be mapped 
out exactly using two-loop perturbation theory . Fur- 
thermore, the bare gauge coupling g 2 is a dimensionful 
quantity, which sets the scale in three-dimensional gauge 
theories. 

We work with the lattice action defined by 

S[U,<j>] =/3 G £;(l-§Trt/ p ) (1) 
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where U^x) G SU(2) is the link variable, U p denotes 
the plaquette, and <fi(x) is the scalar field. The bare 
parameters (3q = 4/ '(ag 2 ), (3h and /3r are the inverse 
gauge coupling, the scalar hopping parameter, and scalar 
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self-coupling, respectively. In our Monte Carlo simula- 
tions we update the gauge fields using a combination of 
the standard heatbath and over-relaxation algorithms for 
SU(2) ]ll[|l2|| . The scalar fields are updated using the 
algorithm described in E3]. Further details about our 
simulation procedure can be found in ref . (^] . 

The standard operator used to compute the static po- 
tential is the Wilson loop, i.e. the correlation of a string 
of length r over a time interval t. Throughout this work 
we have computed Wilson loops of area r x t for "on- 
axis" orientations in the spatial directions, so that all 
separations are integer multiples of the lattice spacing a: 

Gss(r, t) = (Tv {tT(0, rj) U{r% rj + i3) 

xU^(m,rj + m)U^(0,t3)\\, (2) 

where U(x,y) is a shorthand notation for the straight 
line of links connecting the sites x and y. A non-local, 
gauge-invariant operator describing the bound state of a 
static colour source and a scalar field in three dimensions 
has been studied in j[o|, viz. 

G M (f ) = ( \ Tr { 0t ( t 3) f/t (o, t%) 4>(Q) }> . (3) 

This motivates our definition of an operator which 
projects onto two of these bound states separated by a 
distance r 

G MM (r,t) = (iTr{(l-a 3 )^(t3){/t(o,t3)0(O)} (4) 
x Tr {(1 - a 3 )^(rj) {/(rj, rj + i3) 0(rj + «)}). 

Finally, we consider operators which describe transitions 
between a string and a two-meson state and vice-versa: 

G SM (r,t) = (iTr{0t(i3)C/t( O ,t3)[/(O,rj) 

x E/(rj,rj + t3)0(rj + t3)}) (5) 

G MS (r, t) = ( 3 Tr { ^ ( r j) £/(rj, rj + 13) 

x C/t (i3, rj + i3) C/t (0, i3) 0(0) } ^ (6) 

In order to study the static potential and the mixing 
between gauge string and two-meson states, we construct 
the following matrix correlator 

G(r,t)=( r SM t%)- ( 7 ) 

v 1 \G M s{r,t) G M M(r,t) J v ; 

The operator Gjf(t) serves to extract the energy Em of 
a single bound state of a static source and a scalar field, 
which helps to identify the two-meson state, since one 
naively expects Emm — 2Em- 

In order to enhance the signal of all our operators, 
we have constructed smeared gauge and scalar fields. 
Following the standard algorithm in we constructed 



"fuzzed" link variables of length a. Fuzzed scalar fields 
were constructed using a suitable adaptation of the block- 
ing algorithm used in (| [h]]. In particular, eq. (3) of 
ref. H carries over literally after replacing the blocked 
links by their fuzzed counterparts. We constructed spa- 
tial Wilson lines for three different link fuzzing levels, 
whereas five fuzzing levels were used for scalar fields, so 
that we obtained the correlator G(r, t) as an 8 x 8 matrix. 

A procedure of diagonalising G(r, t) followinga varia- 
tional method has been described in detail in ||Jl5|] . Here 
we only wish to point out that the diagonalisation of 
the 8x8 correlation matrix Gjfc yields eigenvectors $^ 
from which the correlation function of the (approximate) 
eigenstates of the Hamiltonian can be calculated, 

8 

r,(t) = ($i(t)$i(0)) = ai^ikG jk (t), (8) 

The coefficients an- quantify the overlap of each indi- 
vidual correlator Gik with the correlators of the mass 
eigenstates, Ti. The diagonalisation procedure allows the 
computation of the ground state energy, the first few ex- 
citations, and, most crucially, their composition in terms 
of the original operator basis. The coefficients play a 
central role for the interpretation of string breaking as a 
mixing phenomenon. 

Our choice of bare parameters followed closely our ear- 
lier work ]^-[lo[|. In particular, we approach the contin- 
uum limit for a fixed ratio PgPr/Ph = 0.0239, corre- 
sponding to a fixed continuum scalar self-coupling. We 
selected three values for the bare gauge coupling, i.e. 
(3g = 5.0, 7.0, 9.0, and the hopping parameter (3h was 
always chosen such that the model was in the confine- 
ment region of the phase diagram. Lattice sizes were 
chosen such that L/j3a — 5.7. At @g — 5.0 we also 
considered L//3g = 7.2 in order to check for finite-size 
effects. The maximum separation of colour charges was 
''max = L/2 in all simulations. At f3o — 5.0 we used 
30 iterations in the fuzzing algorithm, which was suf- 
ficient to observe saturation in the projection onto the 
ground state. The link/staple mixing ratio was always 
set to two. For the runs at smaller lattice spacing, the 
maximum number of fuzzing iterations was scaled by the 
respective ratio of /?G-values. The correlators G(r, t) and 
GmW were measured after every full update, and 50 
individual measurements were collected in bins for post- 
processing. Typically we performed 1500 measurements 
(30 bins) and 3000 measurements for our most accurate 
determinations. Statistical errors were estimated using a 
jackknifc procedure. 

In Fig. [l] we show the energies of the two lowest eigen- 
states as a function of r/a. The ground state shows the 
familiar linear rise for small separations. For r/a > 15 
however, the linear rise saturates, and a comparison with 
the energy extracted for the single meson state shows that 
E(r) ~ 2Em in this region. The qualitative and quan- 
titative behaviour of the ground state energy at large 
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separations is entirely consistent with the string break- 
ing picture in the sense that the ground state energy is 
bounded by that of the two-meson state. Thus, we take 
rb ~ 15a as a rough estimate at this point in parameter 
space. The energy of the first excited state is constant for 
r < n, and consistent with that of a two-meson state. A 
significant r-dependence in this region is only observed 
at very small separations. There are a number of pos- 
sible explanations, such as interactions between the two 
mesons, mixing with the excited potential, or that the 
separation may be too small for a two-meson state to 
exist. By contrast, for r > r*b, the energies of the first 
excited state rise linearly and appear to be the contin- 
uation of the ground state potential below the breaking 
scale. 

aE(r) 
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FIG. 1. The energies of the ground state and the first ex- 
cited state for (3 G = 7.0, j3 B = 0.3468 on 40 3 . The dashed 
lines indicate the location of twice the energy of the single 
meson state, as extracted from Gm(£). 

The interpretation of the r-dependence of the spec- 
trum as clear evidence for string breaking is corrobo- 
rated by the analysis of the composition of the energy 
eigenstates in terms of the original operator basis. In 
Fig. [2] (a) we plot the maximum projection of the string 
and two-meson operator onto the ground state. At small 
r/a it is obvious that the ground state is dominated by 
the string operator, consistent with the picture of static 
quarks bound by a flux tube, whereas the projection of 
the two-meson operator onto the ground state is always 
significantly lower, but clearly non- vanishing. As men- 
tioned above, the admixture of Gmm to the ground state 
may have physical reasons, but it could also arise through 
overlapping fuzzed links in Gmm, which effectively cre- 
ate a loop of links. Operator artefacts of this kind were 
already encountered in Thus, the most likely inter- 
pretation of the overlaps at small r/a is a combination 
of physical effects and operator artefacts. With increas- 
ing r/a Gmm almost fully projects onto the two-meson 



state. For r ~ r^ however, there are clear indications for 
mixing between the flux tube and the two-meson state, 
with comparable projections of both operator types. For 
r > rb the ground state receives almost exclusively con- 
tributions from the two- meson operator, and it is now 
the first excited state (not shown), whose dominant op- 
erator content comes from the string. We conclude that 
there is clear evidence for the crossing of energy levels 
associated with gauge strings and two-meson operators 
at rb, so that the flux tube is energetically disfavoured 
for r > tv 
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FIG. 2. The maximum projections on the ground state of 
the string (ai t ss, full circles) and the two-meson operator 
(ai,MM) open squares) for several values of the hopping pa- 
rameter at Pa = 5.0 on 26 3 . 

At both (3 G = 5.0, /3h = 0.3520 and p G = 7.0, p H = 
0.3468 we observe a fairly sharp crossover and thus a rela- 
tively narrow mixing region. In order to explore the mix- 
ing phenomenon for smaller scalar masses, we have simu- 
lated a sequence of /3jj-values at pa = 5.0 fixed, which all 
satisfy the constraint Pq Pr/ P 2 H = 0.0239, i.e. correspond 
to the same continuum scalar coupling. For smaller scalar 
mass (larger Ph) one expects r^ to occur at smaller sep- 
arations, since less energy is required to support a two- 
meson state. Also, according to || the width of the mix- 
ing region becomes larger for smaller scalar mass, so that 
the crossover is expected to widen as Ph is increased. 
Our plots of the maximum projections as a function of 
Ph displayed in Fig. |^ (a)-(c) show that both is indeed 
the case. The conclusions of ||, which have been ob- 
tained for SU(N) gauge theories with fermions, are thus 
confirmed qualitatively in our simplified model. 

By comparing the spectrum for different Ph at Pq = 
5.0 fixed, we observe that energies extracted from eigen- 
states with dominant admixtures of gauge strings remain 
largely unaffected by the value of Ph, and only the en- 
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ergies of states with a dominant overlap from two-meson 
operators show a significant dependence on (3h- This is 
in line with our earlier observation that scalar and gauge 
sectors decouple approximately JH^]. The consequence 
of this in the present context is that the string operator 
predominantly projects on the state with unbroken flux 
tube for all spatial separations outside the mixing region. 
Hence, Wilson loops exhibit an area law for correlation 
distances t/a up to which the signal can be followed, 
even if r/a is already beyond rt,. This offers a possible 
explanation for the failure to detect string breaking in 
full QCD simulations using only Wilson loops. Namely, 
since the only significant overlap of a string onto the bro- 
ken flux tube occurs inside the mixing region and van- 
ishes again when r/a is increased further, the saturation 
of the energy can only be observed if the mixing region 
is broader than, say, two lattice spacings. Furthermore, 
the width of the mixing region is sensitive to the mass 
of the matter considered. Therefore, in computations of 
Wilson loops in full QCD with relatively heavy dynami- 
cal quarks, the mixing region may be too narrow. If this 
scenario holds, then the problem can only be overcome by 
employing a variational approach of the type discussed 
in this letter. 

We now consider the scaling behaviour of the breaking 
scale rb- As a definition of rb we use the point where 
the maximum overlaps of string and two-meson operators 
onto the ground state are equal, i.e. 



A = ai : ss — 0,1 mm — 0, 



(9) 



where we have borrowed from the notation in cqs. (S), (Q) 
and (||). An estimate of r^/a is obtained through a local, 
linear interpolation of A to the point where it vanishes. 
This procedure becomes exact in the continuum limit. 
Within our accuracy we assign an error of at most one 
lattice spacing to the estimate of r*b/a. Our results for all 
lattice sizes and parameter values are shown in Table |. 
We note again that r^/a decreases as the scalar mass 
is decreased. Furthermore, by comparing the estimates 
for rb/a on L/a = 26 and 36 at 0a = 5.0 we conclude 
that finite-size effects are practically absent in our sim- 
ulations. In the last column rb is given in units of the 
continuum gauge coupling constant g 2 . Within errors we 
observe scaling of the estimates for r^ g 2 from runs c, e 
and /, whose bare parameters are related through lines of 
constant physics in the continuum. For a controlled con- 
tinuum extrapolation, a further value at larger f3a would 
be desirable, so we choose to quote as our final result 



r b g 2 « 8.5. 



(10) 



This number can be combined with results for the mass 
spectrum in the continuum limit at the same physical 
couplings [|. Taking m s /g 2 = 0.839(15) for the lightest 
scalar bound state and mc/g 2 = 1.60(4) for the lightest 
scalar glueball M, we obtain 



To summarise, our results clearly show that string break- 
ing occurs in non-Abelian gauge theories with matter 
fields. The main evidence is the existence of a mixing 
region, where both the string and two-meson operators 
have a significant overlap onto the ground state. We have 
also shown that Wilson loops are not suitable to compute 
the ground state potential for r > rb. Estimates of rb 
are also of phenomenological interest. For instance, rb 
is expected to be related to the scale parameter of QCD 
through light-quark constituent masses ]lq] . Also, r\> is 
important for the understanding and description of frag- 
mentation processes |l7| . An extension of our analysis to 
QCD is therefore highly desirable. 

TABLE I. Estimates for the breaking scale rb. 



Run 


Pa 


Ph 


L/a 


r b /a 


rbff 2 


a 


5.0 


0.3538 


26 


8.2 ±0.5 


6.56 ±0.38 


b 


5.0 


0.3534 


26 


8.9 ± 1.0 


7.12 ±0.79 


c 


5.0 


0.3520 


26 


9.8 ±0.7 


7.88 ±0.59 


d 


5.0 


0.3520 


36 


9.9 ±0.8 


7.90 ± 0.64 


e 


7.0 


0.3468 


40 


14.8 ± 1.0 


8.47 ±0.57 


f 


9.0 


0.3438 


52 


19.1 ± 1.0 


8.50 ± 0.44 
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